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Abstract 

The Metric Traveling Salesman Problem (TSP) is a classical NP-hard op- 
timization problem. The double-tree shortcutting method for Metric TSP 
yields an exponentially-sized space of TSP tours, each of which approximates 
the optimal solution within at most a factor of 2. We consider the problem 
of finding among these tours the one that gives the closest approximation, 
i.e. the minimum-weight double-tree shortcutting. Previously, we gave an ef- 
ficient algorithm for this problem, and carried out its experimental analysis. 
In this paper, we address the related question of the worst-case approxi- 
mation ratio for the minimum-weight double-tree shortcutting method. In 
particular, we give lower bounds on the approximation ratio in some spe- 
cific metric spaces: the ratio of 2 in the discrete shortest path metric, 1.622 
in the planar Euclidean metric, and 1.666 in the planar Minkowski metric. 
The first of these lower bounds is tight; we conjecture that the other two 
bounds are also tight, and in particular that the minimum- weight double- 
tree method provides a 1.622-approximation for planar Euclidean TSP. 



1. Introduction 

The Metric Travelling Salesman Problem (TSP) is a classical combinato- 
rial optimization problem. We represent a set of n points in a metric space 
by a complete weighted graph on n nodes, where the weight of an edge is 
defined by the distance between the corresponding points. The objective of 
Metric TSP is to find in this graph a minimum-weight Hamiltonian cycle 
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(equivalently, a minimum- weight tour visiting every node at least once) . The 
most common example of Metric TSP is the planar Euclidean TSP, where 
the points lie in the two-dimensional Euclidean plane, and the distances are 
measured according to the Euclidean metric. 

Metric TSP, even restricted to planar Euclidean TSP, is well-known to 
be NP-hard [lO|. Metric TSP is also known to be NP-hard to approximate 
to within a ratio 1.00456, but polynomial-time approximable to within a 
ratio 1.5. Fixed-dimension Euclidean TSP is known to have a PTAS (i.e. 
a family of algorithms with approximation ratio arbitrarily close to 1) [l|; 
this generalises to any metric defined by a fixed-dimension Minkowski vector 
norm. 

Two simple methods, double-tree shortcutting ^] and Christofides' |1, 



131 ]. allow one to approximate the solution of Metric TSP within a factor 



of 2 and 1.5, respectively. Both these methods belong to the class of tour- 
constructing heuristics, i.e. "heuristics that incrementally construct a tour 
and stop as soon as a valid tour is created" 0]. In both methods, we build 
an Eulerian graph on the given point set, select an Euler tour of the graph, 
and then perform shortcutting on this tour by removing repeated nodes, 
until all node repetitions are removed. In general, it is not prescribed which 
one of several occurrences of a particular node to remove. Therefore, the 
methods yield an exponentially-sized space of TSP tours (shortcuttings of a 
specific Euler tour in a specific Eulerian graph), each of which approximates 
the optimal solution within at most a factor of 2 (respectively, 1.5). 

The two methods differ in the way the initial weighted Eulerian graph 
is constructed. Both start by finding the graph's minimum-weight spanning 
tree (MST). The double-tree method then doubles every edge in the MST, 
while the Christofides method adds to the MST a minimum-weight matching 
built on the set of odd-degree nodes. The weight of the resulting Euler tour 
is higher than the optimal TSP tour at most by a factor of 2 (respectively, 
1.5), and the subsequent shortcutting can only decrease the tour weight. 

While any tour obtained by shortcutting of the original Euler tour ap- 
proximates the optimal solution within at most a factor of 2 (respectively, 
1.5), clearly, it is still desirable to find the shortcutting that gives the clos- 
est approximation. Given an Eulerian graph on a set of points, we will 
consider its minimum-weight shortcutting across all shortcuttings of all pos- 
sible Euler tours of the graph. We shall correspondingly speak about the 
minimum-weight double-tree and the minimum-weight Christofides methods. 

Unfortunately, for the general Metric TSP (i.e. an arbitrary complete 
weighted graph with the triangle inequality), the corresponding double-tree 
and Christofides minimum-weight shortcutting problems are both NP-hard. 
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(a) Minimum spanning 
tree 



(b) Depth-first double-tree (c) Absolute minimum- 
tour weight tour 



Figure 1: The depth-first double-tree method: a lower-bound instance 



The minimum-weight double-tree shortcutting problem was also believed for 
a long time to be NP-hard for planar Euclidean TSP, until a polynomial- 
time algorithm was given by Burkard et al. In 0], we gave an improved 
algorithm running in time 0(4'^n^), where d is the maximum node degree 
in the rooted minimum spanning tree (e.g. in the non-degenerate planar 
Euclidean case, d < 4). In contrast, the Christofides version of the problem 
remains NP-hard even for planar Euclidean TSP flT]. 

A natural question about the properties of the two approximation meth- 
ods and their variants is whether the approximation ratios 2 and 1.5 are 
tight, i.e. whether there is a problem instance where the approximate so- 
lution has approximation ratio 2 (respectively, 1.5), or a family of problem 
instances where the approximate solutions approach these ratios arbitrarily 
closely. 

For the minimum-weight double-tree method, the answer to this ques- 
tion is unknown, as observed e.g. in The only existing lower bounds for 
the double-tree method apply to a shortcutting that is performed in some 
suboptimal, easily computable order. An example of such an order is depth- 
first tree traversal; we shall call the resulting method depth-first double-tree 
shortcutting. A tight lower bound for this method is given by the stan- 
dard Euclidean lower-bound construction shown in Figured], whichadapted 
from Figure [la] shows an instance point set and the (unique) mini- 
mum spanning tree. We assume that e = o(l); for example, we can take 
e = 1/n. The vertical size of the instance set is 1, and the horizontal size is 
(l -|- o(l))n. The weight of the unique MST is (2 -|- o(l))n; the double-tree 
weight is (4 -|- o(l))n. The double tree undergoes no significant shortcut- 
ting, and the resulting tour (Figure [Tb]) still has weight (4 -|- o(l))n. The 
absolute minimum- weight tour (Figure [Tc]) has weight (2-|-o(l))n, therefore 
the approximation ratio on the given instance set is 2. 

For the minimum-weight Christofides algorithm, a tight lower bound is 
given by the standard Euclidean lower-bound construction shown in Fig- 
ure [2] which is adapted from fb] and uses the same conventions as Figure [1] 
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(a) Minimum spanning (b) Minimum-weight (c) Absolute minimum- 

tree Christofides tour weight tour 

Figure 2: The minimum-weight Christofides method: a lower-bound in- 
stance 



The minimum spanning tree has exactly two odd-degree nodes, therefore 
the additional matching consists of a single edge. The resulting Eulerian 
graph (Figure I2bp is already a Hamiltonian cycle, hence no shortcutting 
is required. The weight of the cycle is (3 -|- o{l))n. As before, the abso- 
lute minimum-weight tour (Figure [2c|) has weight (2-|-o(l))n, therefore the 
approximation ratio on the given instance sets is 1.5. 

In the rest of this paper, we address the question of the worst-case ap- 
proximation ratio for the minimum-weight double-tree shortcutting method 
in some specific metric spaces. In particular, we give a lower bound on 
the approximation ratio in the discrete shortest path metric^; this bound 
is tight, and can be regarded as a lower bound for a generic metric space. 
We also give the first non-trivial lower bound for the planar Euclidean and 
planar Minkowski metrics. 

2. The lower bounds 

2.1. The discrete shortest path metric 

The worst-case approximation ratio of the double-tree and Christofides 
methods can clearly be dependent on the metric in which the TSP problem 
is defined. In the Introduction, we described a tight lower bound of 2 on the 
worst-case approximation ratio in the planar Euclidean metric, both for the 
depth-first version of the double-tree method and for the minimum-weight 
Christofides method. In contrast, no non-trivial lower bounds have been 
known, to our knowledge, for the minimum-weight double-tree method in 
any metric. A tight bound in the Euclidean metric seems difficult to obtain; 
however, it can be established that the upper bound of 2 is tight in some 
non-Euclidean metrics, and therefore is tight for the generic Metric TSP. 



^The same resuh has been obtained independently by Bilo at al. [^. 
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(a) Minimum spanning (b) Minimum-weight (c) Absolute minimum- 

tree double-tree tour weight tour 



Figure 3: The minimum-weight double-tree method: a lower-bound instance 
in the discrete shortest path metric 

Given a weighted undirected graph, consider the discrete shortest path 
metric on its node set. The distance between two nodes in this metric is 
defined as the weight of the shortest path connecting them in the graph. Let 
n be a power of 2. Let T„ be a rooted tree on n nodes, where the root has a 
single child, which branches off into a complete binary tree with n/2 leaves. 
We construct an instance graph on 2n nodes as follows. First, we create two 
copies of the tree keeping track of corresponding pairs of nodes (i.e. pairs 
of nodes which are copies of the same node in r„). We then give all the edges 
in each tree weight 1, and connect the roots of the two trees by a root edge 
of weight 1. Finally, we connect every pair of corresponding non-root nodes 
in both trees by a cross-edge of weight 1 -|- e. We assume that e = o(l); for 
example, we can take e = 1/n. The instance graph corresponding to n = 8 
is shown in Figure [Sal where edges of weight 1 and 1 + e are represented, 
respectively, by solid and dotted lines. 

The unique MST consists of both copies of T„ plus the root edge, and 
has weight (2 — o(l))n; the double-tree weight is (4 — o(l))n. Note that for 
any two nodes o, b within the same copy of the distance between a and b 
is equal to the weight of the path connecting these nodes in the tree. Hence, 
a shortcutting from a, b, c to a,c can reduce the tour weight, only if a and c 
belong to different copies of T„. Also note that any double-tree Euler tour 
of Tn has weight 2n — 2. Any Hamiltonian cycle of the complete weighted 
graph obtained by shortcutting the double-tree Euler tour will contain a 
Hamiltonian path in a complete weighted subgraph induced by each copy 
of T„. The weight of a such a Hamiltonian path can differ from the weight 
of the double-tree tour of T„ by at most the weight of a single edge, which 
cannot exceed 21ogn = o{n). Therefore, the resulting Hamiltonian cycle 
still has weight (4 — o(l))n. 
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(a) Minimum spanning (b) Minimum-weight (c) Absolute minimum- 

tree double-tree tour weight tour 



Figure 4: The minimum-weight double-tree method: a lower-bound instance 
in the Euclidean and Minkowski metrics 

The minimum-weight double-tree tour for our example is shown in Fig- 
ure [SbJ where straight edges have weights 1 and 1 + e, and curved edges 
have integer weights greater than 1. An edge's curvature indicates the lay- 
out of the shortest path along which the edge weight is measured. The 
absolute minimum- weight tour has weight (2 -|- o(l))n, and consists of the 
root edge and all the cross-edges, linked together by edge-disjoint paths in 
the two trees. The absolute minimum-weight tour for our example is shown 
in Figure [3cl using the same graphic conventions as in Figure I3b[ The ap- 
proximation ratio of the minimum-weight double-tree method on the given 
instance set is 4/2 = 2, which matches the generic upper bouncd. 

2.2. Euclidean and Minkowski metrics 

Compared with the above construction for the discrete shortest path 
metric, it appears to be much more difficult to obtain a tight bound in planar 
Euclidean- type metrics. We describe a construction that provides the first 
non-trivial lower bound on the approximation ratio of the minimum-weight 
double-tree method in the planar Euclidean and Minkowski metrics. 

The proposed construction consists of 6n + 1 points, and is shown in 
Figure Hal for n = 4. The instance point set consists of seven points forming 
a symmetric three-way central star of arbitrary constant size, and six rows of 
points extending from the star's ends in three symmetric directions in steps 
of length 1. Figure Ha] shows the (unique) minimum spanning tree, which 
has weig ht [<6 + o{l))n. Fi gure I4bl shows the minimum- weight double-tree 



^Many variations on the described construction are possible. We have chosen a variant 
that is easy to visualise. 
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tour, which has weight (8 + -^3 + o(l))n. Figure He] shows the absolute 
minimum- weight tour, which has weight (6 + o(l))n. The approximation 
ratio of the minimum-weight double-tree method on the given instance set 
is (8 + \/3)/6 « 1.622. There remains a substantial gap between this lower 
bound and the generic upper bound of 2, which is also the best known upper 
bound in the planar Euclidean metric. 

The same construction provides a somewhat stronger lower bound in a 
metric defined by the hexagonal norm — a Minkowski vector norm with 
the unit disc in the shape of a regular hexagon (see Figure Hal) . In this 
metric, the distance between two points is measured along a polygonal path 
composed from segments parallel to the edges of the unit disc. The weights 
of the minimum spanning tree (Figure Ha|) and of the absolute minimum- 
weight tour (Figure l4cll on the above instance set remain asymptotically 
unchanged in the new metric. However, the weight of the minimum-weight 
double-tree tour (Figure l4b|) increases to (lO + o(l))n. Therefore, the lower 
bound in the hexagonal metric is 10/6 ~ 1.666. 

3. Conclusions 

In the previous section, we presented lower bounds on the minimum- 
weight double-tree method. We have shown that the trivial upper bound of 
2 is tight in at least some metrics (in particular, the discrete shortest path 
metric). However, in the important cases of the Euclidean and Minkowski 
metrics, a substantial gap remains between our lower bounds of 1.622 (re- 
spectively, 1.666) and the trivial upper bound of 2. Considering the appar- 
ent difficulty of improving on these lower bounds, and the good approxi- 
mation behaviour of the minimum-weight double-tree algorithm on typical 
Euclidean TSP instances we conjecture that these lower bounds are 
tight, and that the minimum- weight double-tree method provides a 1.622- 
approximation for planar Euclidean TSP. 
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